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Abstract 

Using an approach based on many body perturbation theory, the correlation energy f co is ex- 
pressed as an explicit functional of p±, v, and v s , where p\ is the one-particle density matrix from the 
noninteracting, or reference, determinantal-state; v is the external potential from the interacting, 
or target, state; v s is the (kernel of the) external potential from the noninteracting determinantal- 
state. In other words we have S co [pi,v,v s ]. Anther possibility is the following explicit functional: 
£co [pi i v co j v s ] , where v co is the (kernel of the) correlation potential from the noninteracting Hamil- 
tonian. The proposed method can, in principle, be used to compute £ co in a very accurate and 
efficient manner, since, like the Kohn-Sham approach, there are no virtual orbitals to consider. 
However, in contrast to the Kohn-Sham approach, £ co is a known, explicit functional that can be 
approximated in a systematic manner. For simplicity, we only consider noninteracting closed-shell 
states and target states that are nondegenerate, singlet ground-states; so, in that case, p\ denotes 
the spin-less one-particle density matrix from the determinantal reference state. 
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I. INTRODUCTION 



The Kohn-Sham version of density functional theory has been very successful in 
the description of electronic structure for quantum chemistry and condensed matter 
physics.— i 2 i 3 i 4 i 5 i 6 i 7 Unlike pure density functional approaches, the Kohn-Sham method 
uses a set of occupied orbitals from a noninteracting state, where this determinantal state 
shares its electron density p with the target, or interacting, state In the Kohn-Sham ap- 
proach, the kinetic energy (through the) first order is not an explicit functional of p, but, in- 
stead, this functional depends on the one-particle density matrix pi from the noninteracting, 
determinantal state. The Kohn-Sham approach requires as input the exchange-correlation 
functional E^ — or equivalently the exchange E x and correlation-energy E^ functionals — 
where E^p is required to be an explicit functional of p. Unfortunately, £^ is an unknown 
functional, and there is no systematic method to improve approximations. 

The optimized potential metho d 9 ! 10 ! 11 ! 12 ' 13 ' 14 ' 15 is a density functional approach that can 
convert a nonlocal operator into a local one, where the exchange-correlation functionals can 
depend on both the occupied and virtual orbitals. Unfortunately, this method lacks the 
efficiency of other Kohn-Sham approaches. Furthermore, the optimized potential method 
introduces functional that — in contrast to many wave function methods — are not invariant 
to a unitary transformation of either the occupied or virtual orbitals; the local potentials 
also depend on the orbital energies. 

High levels of approximations — beyond the Kohn-Sham approaches — can be obtained by 
wave function methods , 16 i 17 i 18 i 19 including the coupled cluster method, many-body pertur- 
bation theory, and configuration interaction. Often, however, these methods are much less 
efficient than the Kohn Sham approaches, where, typically, wave function methods consider 
a large number of 2-electron molecular integrals, depending on both the occupied and virtual 
orbitals, and these integrals must be computed and utilized in calculations involving large 
atomic-orbital basis sets. 

In order to improve the efficiency of the wave function methods, especially in regards 
to their scaling with molecular size, perturbative methods have been developed based on 
localized molecular orbitals . 20 ' 21 ' 22 ' 23 ' 24 ' 25 ' 26 An alternative approach uses a Laplace transform 
to remove the energy denominators in perturbation theory, yielding approaches involving 
correlation energy expressions that depend explicitly on the atomic-orbital basis set . 27 ' 28 
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This method has also been used with the coupled cluster theory. 29 

Our research interest is in the development of methods that bridge the gap between 
wave function approaches and Kohn-Sham density functional theory. However, these two 
different approaches already share some common features. For example, as mentioned above, 
Kohn-Sham introduces a determinantal wave function and treats the first-order kinetic 
energy in the same manner as in wave function methods, 1.6. « cLS cL functional of the one- 
body density matrix from a determinant. Furthermore, hybrid density f unct ionals , 30 > 31 ' 32 > 33 
including B3LYP , 30 i 34 introduce a component of the exact exchange — a functional of p\— 
even though these approaches violate the Hohenberg-Kohn theorem^ by using the nonlocal 
exchange operator. 

A simple generalization of the Kohn-Sham functionals involves using the exact exchange 
energy £ x (with its corresponding non-local operator) and, in addition, a correlation-energy 
functional that also depends on p\. Such an approach could probably be derived within a 
variant of the Kohn-Shan scheme. However, it is also reasonable to base such a formalism 
on wave function methods, since, for example, the correlation-energy, say £ co , from wave 
function methods is an implicit functional of p\, and does not implicitly depend on the 
orbitals, even though there is an explicit dependence. (This is easily proven by noting that 
we have £ co = £—E\[px], where £ is the exact electron energy, and E\ is the first order energy 
that is determined by p\.) Furthermore, such a formalism can be based on any reasonable 
orbitals: Hartree-Fock, Brueckner, and natural, but there is probably no advantage to using 
Kohn-Sham orbitals. In addition, however, when considering a wave function approach, £ co 
also has a natural dependence on the external potential v from the interacting or target 
state (This is easily proven from coupled cluster theory by noting the Fock operator 
F appearance in the coupled equations and noting that F depends on v.) Therefore, it is 
reasonable to have £ co depending on both on v and p%. Since v is a known function defining 
the molecular structure, the inclusion of the v should not produce any difficulties. 

Below we introduce an approach based on many body perturbation theory , where 
the correlation energy £ co is expressed as an explicit functional of pi, v, and v s , where p\ is 
the one-particle density matrix from the noninteracting, or reference, state; v is the external 
potential from the interacting state v s is the (kernel of the) external potential from 
the noninteracting state. In other words we have £ co [pi,v,v s ]. Here, the v s explicit depen- 
dence appears because the individual diagrams from the perturbative expansion depend on 
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the zeroth-order Hamiltonian, or, equivalently, on v s , since — together with the number of 
electrons — v s defines the zeroth-order Hamiltonian. However, unless the perturbation ex- 
pansion is truncated at some order, it is easily proven that S co does not implicitly depend 
on v s , since an exact S co does not depend of the zeroth-order Hamiltonian. (Infinite order 
summations, e.g., the coupled cluster method, also have S co not depending on v s .) 

This method presented below can, in principle, be used to compute S co in an accurate and 
efficient manner, since, like the Kohn-Sham approach, there are no vitual orbitals to con- 
sider. However, in contrast to the Kohn-Sham approach, S co is a known, explicit functional 
that can be approximated in a systematic manner. Furthermore, in contrast to the optimized 
potential method , WW 3 ! 14 ! 15 the functionals are invariant to a unitary transformation of 
occupied orbitals. The proposed method, like the efficient wave function methods mentioned 
above , 27 i 28 i 29 removes the energy denominators by using Laplace transforms, as well as other 
techniques, to obtain the correlation energy as an explicit functional of p\, v, and v s . For 
simplicity, we only consider noninteracting closed-shell states and target states that are 
nondegenerate, singlet ground-states; so, in that case, p% denotes the spin-less one-particle 
density matrix from the determinantal reference state. 

Let us also mention that we have recently derived a generalization of the Kohn-Sham 
approach in which the correlation energy S co is assumed to be an explicit functional of v 
and px . 36 ! 37 ! 38 ! 39 In a similar Brueckner-orbital method, called reference-state one-particle 
density matrix theory , 11 S co is also assumed to be an explicit functional of v and pi, 
where for closed-shell systems, the local density approximation (LDA) and the LYP, BLYP, 
and B3LYP functionals was shown to be compatible with the method.— In both of these 
two approach, errors from Coulomb self-interactions do not occur, nor the need to introduce 
functionals defined by a constraint search. By utilizing the method presented here, explicit 
forms of these two correlation-energy functionals can be obtained.— 

II. THE NONINTERACTING STATE 

A determinantal state is a wave function that is an eigenfunction of a noninteracting 
Hamiltonian, say Hnv s , where Hnv s has no electron-electron interaction terms. These non- 
interacting Hamiltonian Hnv 3 are completely specified by the number of electrons N and 
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the kernel of an external operator, say v s . Explicitly we have 

N 

h Vs = J2£M, (i) 

i=l 

where 

K = -\V 2 + v s , (2) 

and the nonlocal operator v s — which we require to be spin- free — is defined by its kernel v s , 
i.e., we have 

VsXi r i) = J dr 2 v s (r h r 2 )x{r 2 ). (3) 
Using second quantization, our noninteracting Hamiltonian, given by Eq. (jTfl. becomes 

K = 52(i\K\j)%i, ( 4 ) 

where the spin-adapted excitation operator is given by 

Eij = ^2 a Laja, (5) 

and the one-electrons integrals are spin- free and are written in chemist's notation;— these 
integrals use a spatial orbital set, say {x}; this set has the following form: 

V»«r(r,w) = Xt(r)o-(w); <r = «,/?, (6) 

where the spatial and spin coordinates, r and u>, can be denoted collectively by x, e.g., 

ip ia (x) = <ip ia (r,u). (7) 

It is well known that there is a one-to-one mapping between determinantal states and 
their one-particle density matricesp^ say 7. For a closed-shell state we have 

7 (xi,x 2 )= Yxx(ri)a(uJi)x* x {r2)(T*M = Y X*'(riM(r 2 )<W, 2 , (8) 

and the spin-less one-particle density-matrix is 

Pi(ri,r 2 ) = ^7(ri,w, r 2 ,u) = 2 x*(ri)x*.( r 2), (9) 

" ^e{xo} 

where the sum is over the occupied orbitals from Eq. (JHJ); this set of orbitals is denoted by 
{Xo}- Comparing Eqs. (JBJ) and yields 

T(xi,x 2 ) = -pi(r 1 ,r 2 )S mu]2 . (10) 



Since our closed-shell determinantal states are determined by pi, henceforth, these kets are 
denoted by |pi). 

Using Eq. (J3J) and the determinantal state |pi), our eigenvalue problem can be written as 

H Va \Pi)=2[ E e»)|Pi> = S>|Pi>, (11) 
Ve{xo} / 

where the set of occupied orbitals {Xo} is defined by the lowest N eigenstates of jF Va , 

•F' v 3 Xw ^wXwi Xw ^ {Xo} ¥ Pi, (12) 

and the right side indicates that these orbitals, from the set {Xo}, determine p\. The 
unoccupied, or virtual, orbital set {Xu} are also obtained from T Vs : 

FvsXr = trXr, Xr £ {Xu} > Pi, (13) 

where these orbitals also determine pi, since, for a complete set of orbitals, the unoccupied 
set determines the occupied set; two sets of orbitals that differ by a unitary transformation 
are considered equivalent. 

While there is no restriction placed upon v s , normally it contains the external potential 
v from the interacting state \^nv), the Coulomb potential J, and an additional term, the 
exchange-correlation potential v xc , that we also permit to depend on the one-particle density 
matrix: 

v a = v + J p . + v^ = v + J Pa + K P1 + vg, (14) 

and the latter expression defines the correlation potential £>£J, where K is the exchange 
operator, defined below; the Coulomb potential J, given by 



J P s(ri) = j r 12 1 p s (r 2 )dr 2 , (15) 
depends on the electron density p s of the noninteracting state |pi): 

p s (r) = pi(r,r) = 2 E Xw(*)x*w(*), Pi — > Ps, (16) 

and we use Eq. @; the right side of this equation indicates that p s is determined by p\. 

The non-local operator v xc can be identified as the exchange-correlation operator. If |pi) is 
the Kohn-Sham noninteracting state, this operator is local and it implicitly depends only on 
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p s , the density from \pi) and the interacting state \^nv)- F° r Brueckner orbital theory , 40 i 45 
this operator is nonlocal and may depend on the external potential v. For Hartree-Fock 
theory, the correlation portion of v xc is absent; so, in that case, v^l is just the exchange 
operator K, where the kernel of this operator is 

K Pl (r 1 ,r 2 ) = -^rr 2 Vi(ri,r 2 ). (17) 

For latter use, we also write down the expression for the Fock operator: 

F Pl = -\V 2 + v + J Ps +K pi . (18) 



III. THE INTERACTING STATE 

Now consider an interacting system of electrons. The electronic energy, say S^ VJ for an 
eigenstate of an interacting Hamiltonian is determined by the number of electrons N and the 
one-body external potential v. This must be the case, since the non-relativistic, electronic 
Hamiltonian is entirely specified by N and v: 

N N N N 

h Nv = -\ + E v ® + 2EE^> ( 19 ) 

i i ijtj j 

where r^ 1 is the electron-electron interaction operator. Furthermore, since each N and v 
yields a different ground state wave function | ^ at^ ) Therefore, the ground state energies 
£ Nv , given by 





Hnv\ 




(®Nv\ 





(20) 



are unique functionals of iV and v. 

Using a second quantization approach, the Hamiltonian does not depend on N, and it 
can be expressed by 

H v = 5>iH v2 )b')4- + X>>|J A + \ *£mi)E ijk u (21) 

ij ij ijkl 

where the symmetry-adapted excitation operators are given by Eq. © and 

Ejjkl = 2j a L a fcA Q ^%o-' ( 22 ) 

and the two electrons integrals are spin-free and are written in chemist's notation.— 
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IV. THE CORRELATION ENERGY AS IN EXPLICIT FUNCTIONAL OF V, p\ 

AND V s 



Using many body perturbation theory , 16 i 17 i 18 i 35 we partition our Hamiltonian into the 
zeroth order part and a perturbation V: 

H v =H Va + V Pl , (23) 

where we choose H Vs , from Eq. as our zeroth-order Hamiltonian. Using Eqs. (j2J, (JIJ), 
(EU, and we have 

^ = \ Y^(ij\kl)E m + Y^(i\v\3)Eij ~ X>k|j)4'- (24) 

ij'fci ij 

The energy through the first order — the zeroth-order energy plus the correction from the 
1st order — for many body perturbation theory is independent of the zeroth-order Hamilto- 
nian. This energy functional, given by 

E 1 \p 1 ,v] = (p 1 \H v \p 1 ), (25) 

depends on the two functions: the external potential v and the one-particle density matrix 
P\. Explicitly, this functional is given by 

E l [p u v) = J dv x [-|v>(r l5 r 2 )] r2=ri (26) 
+ J drx v (ri)p(ri) + ~ J J r^dr 1 dY 2 p{r 1 ) p(r 2 ) - J J J r^dr 1 dr 2 pi(r 1 , r 2 ) pi(r 2 , r^. 

In many-body perturbation theory, the electronic energy is given by an order-by-order 
expansion 

£ = £ + £i + £2 + £3 + ■ ■ ■ , (27) 
where the energy through the first order is 

E 1 = £ + £ u (28) 

and the correlation energy £ co includes everything beyond the first order: 

£ co = £ - E 1 = £ 2 + £ 3 + ■ ■ ■ . (29) 
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In order to obtain the an order-by-order expansion of the electronic energy S, it is con- 
venient to write the perturbation V Pl , given by Eq. (}2"2|) , in normal ordered for m 16 ' 17 using 
|px) as the vacuum state: 

V P1 =V C ^ + Vr+V 2 p \ (30) 
where the constant, one-, and two-body terms are given by 

V? = E 1 = E 1 -E , (31) 
Vf 1 = (32) 

ij 

V 2 P1 = lj2(m){E m } n , (33) 

ijkl 

where the correlation potential is defined by Eq. (JHJ), i.e., 

vg = v.-v-J Pa -k in , (34) 

and £ , Ex, J, and K are given by Eqs. (fTTj) . (J2T?)) . (JT3J), and (JUJ), respectively; the brackets 
{• • • } pi indicate that the excitation operators are in normal ordered form with respect to 
the vacuum state \pi). Furthermore, subtracting Eq. (|T5|) from (0) and using Eq. we 
have 

v^=F pi -F pi . (35) 

Eq. (|3^j) is the most convenient expression for V\. However, we can also substitute Eq. 
into Eq. (}3*2*j) . yielding individual diagrams from the perturbation expansion that depend 
explicitly on v and v s . Substituting Eq. (j33|) into Eq. separates V\ into the zeroth-order 
part T and the Fock-operator part F. (Note that {F} Pl is the one-body portion of the 
Hamiltonian.) 41 While the individual terms from the perturbation expansion can depend on 
v s or F, the correlation energy S co , unless approximated from a truncated expansion, cannot 
depend on either v s or F; wg now prove t licit £co 

is an implicit functional of v and p\. 
Consider the electronic energy S, determined by v and N, while E\ is an explicit functional 
of v and p\. Hence, the correlation energy, given by (£ co = £ — E\), is determined by v, pi, 
and N, or just v and pi, since p\ determines N. Therefore, as mentioned above, at least for 
the ground state, the correlation energy is an implicit functional of v and pi, i.e, £ co [v,pi\. 
We now demonstrate that the individual perturbation terms, or diagrams, can be written 
as explicit functional of v, v s and p\. 
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FIG. 1: The three diagrams from the second-order correlation-energy f co and two other diagrams. 

Generally speaking, perturbation theory expresses S co as a functional of the occupied and 
unoccupied orbitals from \pi), their orbitals energies {e}, and the one-body portion or the 
perturbation V±, which, in our case is — "O^o; ( v ~ v s + J Ps + K Pl ), or (F — JF)), as indicated 
in Eqs. ()32|) . (|34p. and (|35|). The correlation-energy diagrams, however, do not depend on 
the constant portion of the perturbation V c ; also, the two-body part Vi is universal, since 
this operator is the same for all electronic systems. 

Consider for the second-order (correction to the) energy, given by 

£2 = £211 + &2d + &2x-, (36) 

where the Goldstone diagrams for these terms are given by the first three entries of Fig. Q 
respectively. For the first term, we have 

£2uiXo,Xu,e,v£} = 2 E ^^(H^>)HW, (37) 

where this diagram depends explicitly on — the kernel of v{£ — the orbital energies {e}, 
and both the occupied and unoccupied orbitals, {x } and {Xu}- Writing out the integrations 
explicitly, we have 

^2u[Xo,Xu,e,v^} = I I dr 1 dr 2 — ^^ 1 (r 1 )x„,(ri)x^(r 2 )^ 1 (r 2 )xr(r2)x*(ri). (38) 

In order to remove the energy denominator, we use the following Laplace transform : 19 ^ 7 

1 r°° 

' dte- t( - ew ~ er \ (39) 



J 

giving 



£ 2u [Xo,Xu,e,Vco} = (4°) 

p p poo 

2 J2j J dr ^ J Q ^ie tier x;(ri)C(ri)e-* ie -x w (ri)x;(r 2 )C(r 2 ) X ,(r 2 ). 
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Using Eqs. (fT2j) and (fT3j) . and the observation that ^ s depends explicitly on the kernel v s , 
and, according to Eq. (J53j) . £>£J depends explicitly on v and u g , we have 

£2u[Xo,Xu,V,V s ) = (41) 
tor 

where denotes .7\, s (rj) and the round brackets (• • •) serve to remind us that e' 1 ^ acts 
exclusively upon X*( r i) an d e - * 1 ^ acts upon Xw{ r i)- Note that, since v s and — |V 2 do not 
commute, the exponential operator, defined by 

t 2 

e tT vs = exp v 2 + t),)] =l + t (-±V 2 + v s ) + - (-±V 2 + v.) (-±V 2 + v.) + ■ ■ ■ , 

(42) 

cannot separate into two factors^ 

exp [t(-|V 2 + v s )] ^ exp [t(-|V 2 )] X exp [tv s ] , if [-±V 2 , v s ] ^ 0. (43) 

In order to make our formulas less cluttered, especially when considering high order 
diagrams, we introduce a notation where the brackets are understood: 

$*u\Xo,Xu>v,v 9 ]= (44) 

p r poo 

2 E / / dr ^ / ^ie 4l ^x;(ri)C(ri)e-^X-(ri)x;(r2)C(r2)Xr(r 2 ), 

wr ® 

where T Vs is identical to F Vs , except that it acts exclusively on x\ it does not act on the 
correlation operator vQ (or its kernel v{£) or the complex conjugate function x* ■ Similarly, 
JF*^ is identical to !F Vs , except that it acts exclusively on it does not act on v<£ or x- The 
following identities clarify these definitions: 

e -^ x *( ri )D S ( ri ) x ( ri ) = x^ri^ir^e-^xin), (45) 
e-* 1 ^x(ri)uS(ri)x*(ri) = x(ri)«S(ri)e _tl ^X*(ri). (46) 

In order to further compress our notation, we define the following transformed correlation 
operator: 

^ = J&.i%(T i )e* j $. i (47) 
where Zj is a combined index of and tj, i.e., Zj = £j. This notation yields 

^zut/oi^pu^uj = ^ y y rfrx(ir 2 y ^ 1 ^ i 1 pi(ri,r 2 )^ 1 (r 2 )/« pi (r2,r 1 ), 
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where p\ is given by Eq. (JSJ), while its orthogonal complement, K Pl , is given by 

K w (r 1 ,r 2 )=2 53 Xr(ri)x^r 2 ), (49) 

and the sum is over the unoccupied orbitals {Xu}- For a complete set of basis functions, we 
have 

<5(n-r 2 ) =X)xi(ri)x?(r 2 ), * e {x } U (50) 



and therefore 

1 
2 

So we can remove k w from Eq. ()48j1 . giving 



5(rx - r 2 ) = - (pi(rx,r 2 ) + ^(r^ r 2 )) . (51) 



p p poo 

£2u[puv,v s ] = J JdndrzJ dt 1 ^ i 1 p 1 (r 1 ,r 2 )C(r 2 )5(r 2 -r 1 ) (52) 

I J J dr ^J o ^i^op 1 1 Pi(ri,r 2 X 1 (r 2 )pi(r 2 ,r 1 ). 



2 

Hence, we have £ 2u as an explicit functional v, v s , and p x . 

In order to evaluate Eq. (|52|) . one must handle the operator exp(£jF) acting on the Dirac 
delta function 5, since the first term on the right side is given by 

r r poo 

J Jdx x dx^J dhv^lp^^v^r^-v,) (53) 
= J j dr x dr 2 J dM£(ri) (e^^Oi, r 2 )) t>£(r 2 ) (e*^(r 2 - r x )) . 

Using the expansion for e tjr given by Eq. (f4*2*j) . this factor will yield terms containing the 
Laplacian operator V 2 acting on 8, i.e., V 2 8, as well as other higher order terms, e.g, V 4 <5. 
One way to handle these terms is by using a basis set expansion for 8, as in Eq. ([5D|h where 
any complete basis can be used that satisfy the boundary conditions. Another possibility is 
to use one of the integral representations of the delta function, for example, 

<J(r a - r 2 ) = -L / e ik < r ^dk. (54) 

OTT J 

Now consider the second term of Eq. (|36p. given by the second diagram from Fig. 
evaluating this diagram yields the following expression: 

^2d[Xo,X«,e] = 2 53 53 e _ ~ l +e _ e (rw\sx)(wr\xs), (55) 

w£{Xo} x£{Xo} r£{xu} s£{xu} W 
12 



and this term can written as 



£ 2 d[Xo,Xu,e] = 2^2^2(e w - e r + e x - e s ) 1 (56) 

wx rs 

Using a generalization of Eq. (|39|). we have 

£2d{Xo,Xu,t] = 2 ^2^2 J J J J rfr i rfr 2 rfr 3^ r 4 J dt 12 (57) 



xe' 12 ^e' 12e3 rr 2 1 e-* 12e -e-' 12 ^ X ;(ri)x ? ,(ri)x:(r 2 )Xx(r 2 )x;(r3)Xr(r3)r3- 4 1 x;(r4)x s (r4). 
Using our notation, defined by Eqs. and (JUIJ), and Eqs. and Q, we obtain 




£2d[pi,v s ] = - / / I dY^r^Yzdr^ I dt 12 (58) 







x r 



12) 



Pi(ri, r 3 )K pi (r 3 , ri)pi(r 2 , r 4 )/t pi (r 4 , r 2 )r 34 1 , 



where the transformed electron-electron interaction term (r^) 1 depends on the variables 
t\ii ri, r 2 , and the function v s ; this operator has the following form: 

1_ = e ti2Fi: e t 12 ^; r -i e -tufi e -t 12 j* V3 _ (59) 

r n 

The 2d subscript indicates that £ 2 d correspond to a second order "direct" diagram. The 
second-order exchange term £ 2x , from Eq. ([36)1 . is represented by the third diagram of Fig. 
This diagram can be evaluated in a similar way as the direct one £ 2 dj giving 
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£2x[f>i, v s ] = -j^ J J J J dridr 2 dr 3 dr 4 y dt 12 (60) 

(7i2)~Vi( r i, r 3 )K pi (r 3 , r 2 )pi(r 2 , r 4 )/t pi (r 4 , r^r^ 1 . 

Note that the and £ 2x do not explicitly depend on the external potential v. 

In Fig. Q the dummy indices that are used are denoted in the diagrams. For example, 
the first diagram has zi for the first vertex and r 2 for the second index, where these choices 
are indicated by 1 and 2 in the figure. In order to evaluate this diagram and obtain Eq. (pfKj) 
directly, we start at the bottom vertex and go backwards: The bottom vertex gives a factor 
of — "Ccopji the hole line from vertex 1 to 2 gives a factor of |pi(ri,r 2 ); the top vertex gives 
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— £>ro( r 2)j the particle line from vertex 2 to 1 gives |K pi (r 2 ,r 1 ). Furthermore, there is an 
additional factor of two for the loop, since the spin state is conserved on an orientated path.— 
Similarly, the second two diagrams from the figure can be evaluated to give Eqs. (J58|) and 
where in both diagrams, the lower interactions, with dummy indices ri, r 2 and tu, 
yield the operator (r^) -1 . For the second diagam in the figure, the direct one £24, we have 
the following factors: a |, since the diagram is symmetric, a factor of 2 x 2 for the two loops, 
and a factor of|x|x|xi for the two p\ terms and two k terms. Putting this together 
gives |x2x2x~x~x|x| = ~, in agreement with Eq. (J5SJ) . 

Now consider the fourth diagram in Fig. ^ where for convenience we denote this third 
order diagram by S^, since this is a direct diagram where the middle interaction lies on two 
hole lines. Evaluating this diagram gives 

( yr I zs) (wy \ xz) (rw | sx) 

wx rs yz v tf 

where we have two energy denominators: One corresponding to the lowest interaction, given 
by (e w + e x — e r — e s ), and one for the middle interaction, given by (e y + e z — e r — e s ), where 
we denote the occupied orbitals by w, x, y, and z and the unoccupied orbitals by r and s. 
Using the following identity?^ 

1 



chh_o^^^ [yr\zs){wy\xz){rw\sx) 

^d- z 2-^2-^2^ r e _ e +e _ € )( € _ e +e _ e y ^ L > 

c r ~ c 2 c sA c tu c r ~ c ie c s/ 



00 



£ j £ i JO Jt 2 



dt 2 e- t2 ^- £l) \ dt±e- tl£ \ (62) 



the product of the two energy denominator can be written as 
1 



\€y 6 r -\- € z (-T 



(63) 



dt u e- t34(ey - ew+€z - ex) / dt 12 e -*M(e«-er+e.-e.) _ 

Comparing the right side of this expression with the the fourth diagram in the figure, or with 
Eq. 1)61)1 . we see that the two lowest interactions, with the following two matrix elements: 
(rw\sx) and (wy\xz), are directly related to the energy exponential factors: £12 (e r — £ w + 
e s — e x ) and £34 (e^ — e y + e x — e z ) . Therefore, multiplying Eq. by these two matrix 
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elements gives 

(rw\sx)(wy\xz) 

(Ej/ ~\~ £s)\£w ~\~ 




rfrirfr2<ir3(ir4 (64) 



x 




fQQ POO 

/ dt u / dt 12 (e* 12 ^e il2 ^rr 2 1 e-* 12 ^e-* 12 ^x;(ri)x.(r 1 )x:(r 2 )Xx(r 2 )) 

JO JtzA 

x (e' 34e -e i34 ^r 3 - 4 1 e- i34 ^e-* 34e ^;(r3)x y (r3)x;(r4)x 2 (r4)) 

a?ri<ir 2 <ir3(ir 4 / ^34 / dt\% 
Jo Jt M 

x(^2)' 1 X;(ri)x w (ri)x:(r2)Xx(r 2 )(r 3 l)~ 1 x;(r3)x y (r3)x:(r4)x 2 (r 4 ), 
where we used Eq. ()59|) . Comparing Eqs. ()64|1 and ()61|1. and using Eqs. Q and (J4T?|) . gives 



1 f f f f f f />oo 

^3d[Pi'^] = — J J J J J J dr i dr 2dr 3 dr 4 dr 5 dr 6 J dt u J dt 12 (65) 
x(f^)" 1 (f : 3|)- 1 r^ 1 pi(ri, r 3 )pa(r 3 , r 5 )/c pi (r 5 , ri)pi(r 2 , r 4 )/9i(r 4 , r 6 )K pi (r 6 , r 2 ). 



Note that in this diagram the two bottom interactions, and , are transformed 
into (r^) -1 and (^j) -1 , respectively, while the top interaction, r^ 1 , is not transformed. 
Furthermore, the transformed interaction that is next to the top, {r^l)~ l , with the dummy 
index £34, gives an integration limit of to 00, while the one below that, (^i 2 ) _1 , with the 
dummy index t 12 , has an integration limit of £34 to 00. By considering higher-order diagrams, 
and by generalizing the identity given by Eq. (pj2|) . e.g.,— 

j poo poo poo 

= / dt z e -*2( £ 3-«) / dt 2 e~* 2(£2 - £l) / d^e^ 161 , (66) 

^3^2^1 Jo Jt 3 Jt 2 

it is easily shown that the following rules apply to all diagram: 1) A factor of \p\ and is 
given for each hole and particle line, respectively. 2) Interactions given by v^{yi) and r~j are 
transformed into and (r^) -1 , unless they are the last, or uppermost, interaction. 3) The 
spatial dummy indices, r 1; r 2 , ■ • • , are integrated over. 4) A tj dummy index from and a 
tij index from (r^) _1 are integrated over using the following rules to obtain their integration 
limits: i) If the interaction — either v^opl or (^ij) -1 — i s the second from the top — the next to 
the last one — the integration limit is from to 00. ii) For other interactions, the integration 
limit is from tj to 00, if the interaction above it is Vcopl] while, the integration limit is from 
tjk to 00, if the interaction above it is (fj^) -1 . In addition, we also have the usual rules of 
symmetry and sign factors as well as a factor of two for each loopi£ corresponding to a sum 
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over each spin state, a and 0. Also, according to Eq. the one-body interaction has a 
negative sign associated with it, it is — v^ p l- 

As a final example, consider the last diagram appearing in Fig. ^ where we denote this 
arbitrary (arb) fourth order diagram by £ ar b- Using the rules above we get 



1 f f f f f f f roo roc roo 

£ aT h[pi,v,v s ] = -— J J J J J J J dr 1 dr 2 dr 3 dr 4 dr 5 dr (i dr 7 J dt 5 J dt u J dt n (67) 
x {rf^irll^v^r^piiri, r 3 )pi(r 3 , r 6 )/t Pl (r 6 , ri)pi(r 2 , r 4 )pi(r 4 , r 7 )n Pl {r 7 , r 5 )K pl (r 5 , r 2 ). 



Hence, using the above rules for all diagram, and Eqs. and (Jo"Tj) — for and 



pi 



yields the correlation energy £ co as an explicit functional pi, v, and v s . i.e, £ co [pi, v, v s ]. The 
other two choices are S co as an explicit functional pi, and v s , or, using Eq. (JHEJ), ^ co as 
an explicit functional pi, u, and J 7 pi . i.e, £ co [pi, ^J, where the v dependence arises from 
the Fock operator F Pl , using Eq. (JT5j) . 

It is well known the the correlation energy cannot depend on the zeroth order Hamiltonian 
T py . Therefore, while the expansion £co[pi, v, J- pi ] has individual diagrams that explicitly 
depend of T px , the entire series does not. Certain infinite-order partial summations also yield 
results that are independent of a zeroth order Hamiltonian. For example, if all diagrams 
are included that correspond to a coupled cluster calculation, the final approximation of £ co 
does not depend on the zeroth order Hamiltonian jF pi , that is, we have the following implicit 
functional: £ co [pi , v\ . 
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